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We describe a C++ program that we have written and made available for calculating the evolution 
of interacting scalar fields in an expanding universe. The program is particularly useful for the 
study of reheating and t hermalization after inflation. The program and its full documentation are 
available on the Web at littp://physics. stanford.edu/gfelder/latticeeasy. In this paper we provide 
a brief overview of what the program does and what it is useful for. 
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I. INTRODUCTION 

Studying the early universe requires describing the 
evolution of interacting fields in a dense, high-energy en- 
vironment. The study of reheating after inflation and the 
subsequent thermalization of the fields produced in this 
process typically involves non-perturbative interactions 
of fields with exponentially large occupations numbers 
in states far from thermal equilibrium. Various approxi- 
mation methods have been apphcd to these calculations, 
including linearized analysis and the Hartree approxi- 
mation. These methods fail, however, as soon as the 
field fiuctuations become large enough that they can no 
longer be considered small perturbations. In such a sit- 
uation linear analysis no longer makes sense and the 
Hartree approximation neglects important rescattering 
terms. What we have learned in the last several years is 
that in many models of inflation preheating can amplify 
fluctuations to these large scales within a few oscillations 
of the inflaton field. Moreover, such large amplification 
appears to be a generic feature, arising via parametric 
resonance in single-field inflationary models and tachy- 
onic instabilities in hybrid models. 

The only way to fully treat the nonlinear dynamics 
of these systems is through lattice simulations. These 
simulations directly solve the classical equations of mo- 
tion for the fields. Although this approach involves the 
approximation of neglecting quantum effects, these ef- 
fects are exponentially small once preheating begins. So 
in any inflationary model in which preheating can occur 
lattice simulations provide the most accurate means of 
studying post-inflationary dynamics. 

Over the past several years we have developed a 
C++ program for doing such lattice simulations. A 
number of studies have already been published by 
us and our collaborators using this program. We 
are now making it available on the World Wide 
Web in the hopes that it will stimulate and aid 



further research in early universe cosmology. We 
call the program LATTICEEASY; its websit e is at 
http: / /physics. stanford.edu/gfelder/latticeeasy/ . The 
website for the program has documentation, including 
derivations of all the equations used in the program. 
Here we present a short summary of what the program 
does and what it can be used for. For more details we 
invite you to visit the website. 

Section ^ of this paper gives an overview of what the 
program is and how it works. Section |l| describes the 
evolution equations solved by the program as well as 
the setting of initial conditions. The references section 
gives a list of papers currently published using LAT- 
TICEEASY results. 



II. OVERVIEW 

Although we designed LATTICEEASY for our studies 
of reheating after inflation the program can more gener- 
ally solve the classical equations of motion for interacting 
scalar fields, with or without the effects of the expansion 
of the universe. 

The study of nonlinear scalar field dynamics has a vari- 
ety of applications. We have used it to study parametric 
resonance , the formation of gravitational waves 

, phase transitions and the formation of topological de- 
fects thermalization after reheating 10|, and 
the formation and evolution of disordered chiral conden- 
sates Q . The program could also be used to study the 
classical limit of general nonequilibrium quantum field 
theories. 

Each particular scalar field potential that the program 
solves is encoded in a model file, which is actually a 
header file read in by LATTICEEASY. For example, we 
have created a model file called twof 1dm. h that contains 
all the necessary equations for running the potential 
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^ = ^mV + ^5W. (1) 

These equations include the potential itself and its first 
and second derivatives, all of which have to be provided 
in the model file. The process of creating a model file for 
a particular potential is described in the documentation. 

Aside from the model file, the only other LAT- 
TICEEASY file that the user needs to modify is 
parameters .h, which contains all the parameters for a 
given run of the program. These include the number of 
grid points, the time step, and a number of other general 
variables specific to each run. 

One of these variables controls the expansion used in 
the run. The run can be set to not include expansion, to 
use a fixed power-law expansion, or to self-consistently 
solve the Friedmann equations using the fields in the 
simulation. (This latter option makes sense if the fields 
are assumed to be the dominant energy component in 
the universe, as is typically the case for simulations of 
reheating.) 

LATTICEEASY has built-in routines for outputting 
means and variances of the fields, field spectra, the com- 
ponents of the energy density, overall energy conserva- 
tion, histograms of the fields, and two-dimensional his- 
tograms of pairs of fields. If expansion is being calculated 
then the scale factor and its derivatives may also be out- 
put. Finally there is an option to call additional output 
routines from the model file so that when you create a 
new model you can design output specific to it. The 
file parameters .h includes a list of parameters control- 
ling which of these outputs will be generated as well as 
specifications such as how often to generate output, how 
many bins to use in the histograms, and so on. 

The figures below illustrate some of the output pro- 
duced by LATTICEEASY. The program's output is in 
the form of ASCII files and these plots were generated 
with Mathematica using those files as input. The LAT- 
TICEEASY website includes a set of Mathematica note- 
books for plotting all the output of the program, but this 
can easily be done using any standard plotting software. 
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FIG. 1. The power spectrum of the field x in the model 
V = ■|A(^'* -I- \g^4>^X^- The spectra rise over time as fiuc- 
tuations of the field are produced. Initially the production 
occurs in resonant peaks but these are quickly smoothed out. 




FIG. 2. Two-dimensional histogram of the fields an 
and Gi in the SUSY F-Term hybrid infiation model 
V = A\\(t)\^ (lap -h lap) +4Ajacr-«2|2. The complex field a 
has fallen down to the minimum of its potential at |a| = v. 

Once all of these parameters have been set you simply 
compile and run LATTICEEASY. The code is designed 
to be platform independent and should work with any 
C-l— I- compiler .0 

III. EQUATIONS 

LATTICEEASY uses a staggered leapfrog algorithm 
with a fixed time step. This means that at each step the 
field values / and their derivatives / are stored at two 
different times t and t + dt/2 respectively. The deriva- 
tives are used to advance the field values by a full step dt 
and then the field values are used to calculate the second 
derivatives /, which are in turn used to advance the field 
derivatives by dt. This evolution is done in place, mean- 
ing the newly calculated field values and/or derivatives 
overwrite the old ones. 

This method is stable for second-order differential 
equations provided they have no first derivative terms, 
i.e. / = /(/). As we will see below the evolution equa- 
tions for scalars in an expanding universe do in gen- 
eral contain first derivative terms, so the program uses 
rescaled variables in which these terms are eliminated. 
These rescalings are described briefly below and in detail 
in the documentation. 



*It's worth noting as a warning that the G libraries used 
by gcc in RedHat Linux 5.0 have a bug that occasionally 
causes the program to crash. This is a compiler bug that was 
fixed in subsequent versions of gcc, e.g. the distribution that 
accompanies RedHat 6.0. Moreover this bug has never in our 
experience caused the program to produce incorrect output. 
When it occurs it causes the program to crash during the 
setting of initial conditions. 
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In this section we always use / to denote any scalar 
field and we use units where Mp ((w 1.22 x lO^^GeV) = 
1. 



A. Evolution Equations 

The evolution eqiiations solved by LATTICEEASY 
are simply the Euler-Lagrange equations for scalar fields 
in an expanding universe 



point. The initial field values are given by quantum fluc- 
tuations with a dispersion characterized by 



• dV 
f + 3Hf+ — = 



(2) 



where / is a scalar field and H = ^ in the Hubble param- 
eter. In addition, if expansion is being calculated self- 
consistently the program solves the Friedmann equations 



47ra , „ , 
a=— 3-(p + 3p) 
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(3) 



(4) 



Either of these Friedmann equations would suffice to 
solve for the expansion. (Either Friedmann equation plus 
the field evolution equations implies the other one.) In 
practice LATTICEEASY uses a combination of the two 
chosen for computational convenience. 

The field and spacetime variables used by the pro- 
gram are rescaled from their bare, physical values. These 
rescalings accomplish several things: They eliminate the 
/ term from the field evolution equations, thus making 
the staggered leapfrog method stable. They can simplify 
the equations by eliminating coupling constants or other 
parameters. They can set the field and time values to 
the natural scales of the problem, e.g. by having the 
time variable correspond to the number of infiaton oscil- 
lations. The program has algorithms for determining the 
most convenient rescalings for a given potential. If you 
wish, however, you can set the rescalings manually to 
whatever scales you find most convenient. The program 
does impose, however, a certain relationship between the 
field and spacetime rescalings that ensures that / will be 
eliminated from the equations of motion. See the docu- 
mentation for details. 



B. Initial Conditions 

Although the field equations are solved in configura- 
tion space with each lattice point representing a posi- 
tion in space, the initial conditions are set in momentum 
space and then Fourier transformed to give the initial 
values of the fields and their derivatives at each grid 



<\M'>=^- 
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where fk is the Fourier transform of / and 
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(6) 



(7) 



The phase of each mode is set randomly and the am- 
plitude of each mode comes from a Gaussian random 
distribution. One of the run parameters is a random 
number seed that can be adjusted to do runs with dif- 
ferent randomized initial conditions. 

In addition to these fluctuations the program allows 
you to set homogeneous initial values for the fields and 
derivatives. The initial field value at each point is a sum 
of the homogeneous value and the fluctuations at that 
point. 



IV. CONCLUSIONS 

There are several additions to the program that we 
hope to make in the future. A great deal of attention 
has been paid recently to the growth of metric perturba- 
tions in the early universe, and we would like to extend 
LATTICEEASY to calculate the coupled equations for 
the growth of field and metric fiuctuations. We also an- 
ticipate adding vector fields to the lattice calculations, 
thus allowing the program to calculate gauge models. 
We have thought about incorporating fermionic fields as 
well, but this presents some difficulties because fermionic 
fields have no well-defined classical limit. At the moment 
it's uncertain if we will find a way to include them in our 
calculations or not. 

The list of references below shows the work that 
has already been published using the results of LAT- 
TICEEASY. As computing power increases and cosmol- 
ogy turns to increasingly difficult mathematical prob- 
lems, the role of large scale numerical simulations in the 
field is likely to keep increasing dramatically. Given our 
uncertainty about the correct particle physics models to 
use for describing the early universe, we believe that the 
most important tools for numerical computation will be 
ones that are flexible enough to easily accomodate a wide 
variety of models. We designed LATTICEEASY to have 
this flexibility while still remaining easy to use and un- 
derstand. 
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